Signatures of Classical Diffusion in Quantum Fluctuations of 2D Chaotic Systems 



Tsampikos Kottos, Alexander Ossipov and Theo Geisel 
Max- Planck- Institut fur Strdmungsforschung und Fakultat Physik der Universitdt Gottingen, 
Bunsenstrafie 10, D-37073 Gottingen, Germany 
(Dated: February 2, 2008) 

We consider a two-dimensional (2D) generalization of the standard kicked-rotor (KR) and show 
that it is an excellent model for the study of universal features of 2D quantum systems with under- 
lying diffusive classical dynamics. First we analyze the distribution of wavefunction intensities and 
$H , compare them with the predictions derived in the framework of diffusive disordered samples. Next, 

we turn the closed system into an open one by constructing a scattering matrix. The distribution of 
the resonance widths V(T) and Wigner delay times V(tw~) are investigated. The forms of these dis- 
tributions are obtained for different symmetry classes and the traces of classical diffusive dynamics 
are identified. Our theoretical arguments are supported by extensive numerical calculations. 



< 



-a 
c 
o 



(N 
> 



PACS numbers: 05.45.Mt, 73.23.-b 



I. INTRODUCTION 

Random Matrix Theory (RMT) was invented more than 50 years ago by Wigner in order to describe the statistical 
properties of complex nuclei In recent years the fast developments in mesoscopic physics and quantum chaos 
gave a new boost to the RMT approach 0, • The strength of RMT consists in the universality of its predictions 
containing no energy or length scales or any parameter dependence. At the same time this is its weakness, because it 
does not allow to take into account various phenomena appearing in mesoscopic systems which introduce new length 
scales or parameters in the system. 

The breakdown of universality and the traces of the underlying classical dynamics are important topics in quantum 
chaos studies and have major applications in mesoscopic physics. It is therefore highly important to compare the 
£j ' results from RMT models, to the results from systems that have a semiclassical limit, and to look for circumstances 
where semiclassical tools rather than RMT become relevant. In this paper we extend our previous studies 0,H| on 
two-dimensional chaotic systems with underlying classical diffusion paying attention to applications of our results for 
random lasers and microwave absorption. We identify the signatures of classical diffusion in quantum fluctuations by 
analyzing the distribution of wavefunction intensities of closed systems, and the resonance widths and Wigner delay 
^— H ' times of their corresponding scattering analogs. 

t-H , Our analysis is based on a two-dimensional (2D) generalization of the well known Kicked Rotor (KR) model . 
For large values of the kicking strength the model is classically chaotic and the dynamics is diffusive. The detailed 
analysis of the classical system allows us to extract the relevant parameters that penetrate into the quantum mechanical 
description. The quantum analog shows dynamical localization, [E 0> H El with eigenstates which are exponentially 
localized in a similar fashion as the eigenstates of disordered systems. However here the localization occurs due to 
complicated interference effects created by the underlying classical chaotic dynamics. Attempts were made to put this 
analogy between the kicked rotor and disordered systems 0, E3 on a solid base but the problem is still not solved 
. completely ||. 

Since the pioneering work of Anderson 11] it is known that the diffusion in disordered systems can be affected very 
strongly by quantum localization. If the system size exceeds the localization length, then the diffusion is completely 
^ ■ suppressed after the time needed for the wave packet to spread over the scale of the localization len gth. But even 
in the absence of strong localization the existence of pre-localized states 0, 0, 0, 0, 0, 0, 0, ITU l20j influences 
the behavior of different physical quantities in the diffusive regime. Having this in mind, it is legitimate to ask, 
' what happens with deterministic chaotic systems when quantum mechanical description becomes relevant? Our 
analysis shows that pre-localized states are present also in the case of dynamical systems and affect various statistical 
quantities like resonances and Wigner delay times. For the analysis of the latter quantities, we turned the closed 
KR to a scattering system by imposing 'absorption' to a set of momentum states. The resulting scattering matrix is 
written in terms of the unitary evolution operator of the corresponding closed system and a projection operator. A 
detailed description of this construction is presented. 

The paper is structured in the following way. In section II the classical model is introduced and the main definitions 
are given. Sections III and IV are devoted to the corresponding quantum mechanical system. The quantization on 
a torus is presented in section III, while in IV we present a detailed construction of the scattering matrix. The next 
two sections V and VI are dedicated to the analysis of various statistical properties of the closed and the open 2D 
KR model correspondingly. Finally, our conclusions are summarized in section VII. 
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FIG. 1: Poincare section of the classical phase space for Hamiltonian Q for (a) k — 0.36, (b) k — 6.37 

II. THE CLASSICAL 2D KICKED ROTOR 

A prominent example of a system with classical diffusion is, the well known in the field of quantum chaos, Kicked 
Rotor (KR) model ;6j, which consists of the free propagating rotor driven periodically in time by an external force. 
Since the total energy of a driven system is not a constant of motion any more, the chaotic behavior can appear even 
in the one-dimensional (ID) case. 

Although many results are known for the standard ID KR there is almost no study for its two dimensional gener- 
alization, besides ref. |3| where the authors have focused on the analysis of dynamical evolution. The corresponding 
model is described by the following Hamiltonian: 

H = H^ + kV{{9 l })^8{t-mT), (1) 

m 

2 



i=l 



where the index i — 1(2) is related to the first (second) rotor correspondingly. Ci denotes the angular momentum and 
9i the conjugate angle of one rotor. The kick period is T, k is the kicking strength, while r is a constant inversely 
proportional to the moment of inertia of the rotor. The parameter 7 is an irrational number whose meaning will be 
explained below. The Hamiltonian Q describes a system which is kicked periodically in time. Another representation 
of the Hamiltonian Q may be given by one rotor moving on a two-dimensional sphere. The two rotors are interacting 
with each other by the potential 

V({9i}) = fcos(0i) cos((9 2 ) cos(a) + i sin(26»i) cos(26» 2 ) sin(a)^ (2) 

The parameter a breaks time reversal symmetry (TRS). The form (|2J was chosen so as to minimize the effect of 
symmetry breaking on the classical diffusion constant D (see bellow). 

The classical equations of motion which correspond to the Hamiltonian Q can be integrated over the kick period 
T giving a set of simple maps: 

0i(n + l) = di(n) +r 1 T£ 1 (n) mod 2tt 
9 2 (n + l) = 9 2 (n) +T 2 T£ 2 (n) mod 2tt 

£1(71 + 1) = d(n) - k(- sin(0i(ra + l))cos(0 2 (« + l))cos(a) + 

+ cos(26»i(n + 1)) cos(26» 2 (n + 1)) sin(a)) 
£3(71+1) = £ 2 (n)-fc(-cos(0 1 (n+l))sm(02(7i + l))cos(a<)- 

- sin(20i(n + 1)) sin(20 2 (n + 1)) sin(a)) (3) 

where 0j(n) and £^(n) are the values of the dynamical variables taken just after the n-th kick. The motion generated 
by this set of maps is integrable in the absence of the kicking potential. For sufficiently small but non-zero k the 
phase space of this system contains both regular islands and chaotic sea (see Fig[TJi). When k is large enough then 
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FIG. 2: Diffusion in momentum space for k = 9.1 

the dynamics becomes fully chaotic (see Fig. and there is diffusion in momentum space (Fig. |2J) with diffusion 
coefficient 



D — limt- 



< L 2 (t) > 



t 



2T 



(4) 



The last expression is correct within the random phase approximation |y, |8|, |21j (see Appendix II) . 

The corresponding open system is described by the same set of equations Q as the closed one, on top of which we 
add absorbing boundary conditions. Namely, we set Li equal to zero, if Li < or Li > L. These conditions give to 
"the particles" the possibility to escape from the system. The evolution of the classical density p(x,y,t) is described 
by the Frobenius-Perron equation, which takes in our case the form of the following diffusion equation: 

dp D 
~dt 

with absorbing boundary condition 



-Ap 



(5) 



Jn — 



"4-(Vp)r 



(6) 



which sets the flux at the boundary J n to be equal to the number of particles <& n carried out under one iteration of 
the map. The solution of Eq. JSJ and Eq. J§J can be represented as the superposition of the diffusive eigenmodes 
Vk{x,y): 



p{x,y,t) 



E 

fe=i 



Cfce lkt v k {x,y) 



(7) 



where 7& are the corresponding eigenvalues and Ck are coefficients determined by the initial condition. The asymptotic 
behavior of the density is governed by the smallest eigenvalue "fi = r c i. As a consequence one has an exponential 
decay of the classical survival probability P(t) — J J dx dy p(x,y,t): 



P{t) oc e~ r ' lt , 



(8) 



The classical decay rate r c ; ~ D /L 2 can be estimated as the inverse time needed for the particle to reach the boundary 
(Thouless time). The exact value of the classical decay rate can be obtained as the solution of the corresponding 
eigenvalue problem or from the numerical calculation of P(t). In Fig. [3] we present the results of our numerical 
calculations for P(t) for some representative parameters. 



III. QUANTIZATION ON A TORUS 

Because of the periodicity of the external force a time dependent solution of the KR model can be represented as 
a superposition of Floquet states [2^J . These are the eigenstates of the evolution operator for one period (Floquet 
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FIG. 3: Exponential decay of the classical survival probability for open kicked rotor model. 

operator). For kicked systems the interaction with external force is instantaneous and one can factorize the total 
Floquet operator into the product of the evolution operators corresponding to the free propagation and the interaction. 
Due to this fact the kicked systems are very convenient for numerical study. 

If the Li are taken mod(27rwi/TTi) where rrii are integers, Eq. JQ) defines a dynamical system on a torus. The 
quantum mechanics of this system is described by a finite-dimensional time evolution operator for one period 

U = exp [-iH Q ({Ci})T/2] exp [-iV({0i})] exp [-iH ({Ci})T/2] (9) 

where we put h = 1. In the above representation, U is symmetric and describes the evolution in the middle of the 
rotations between two successive kicks. Upon quantization, additional symmetries associated with the discreteness 
of the momentum show up, which can be destroyed by introducing irrational values for the parameters 7$'s. The 
most striking consequence of quantization is the suppression of classical diffusion in momentum space due to quantum 
dynamical localization |y,|2|. We introduce the eigenstate components 1 4 r fe(n) of the Floquet operator in the momentum 
representation by 

Y J U m n^k{n)=e iu " T ^ k {n) . (10) 

m 

The quantities tou are known as quasi-energies, and their density is p = T/2-k. The corresponding mean quasi-energy 
spacing is A = l/(pL d ), where L is the linear size of the system and d is the dimensionality (d = 2 in our case). The 
Heisenberg time is tn = 2n/A while to = L 2 /D is the diffusion time (Thouless time). Now one can formally define 
a dimensionless conductance as 

9 = t H /t D = D k L d - 2 (11) 

where D k — TD is the diffusion coefficient measured in number of kicks. Four length scales are important here: the 
wavelength A, the mean free path Im, the linear extent of the system L, and the localization length £. In the rest of 
this paper we will always assume that 



A < l M < L < £. 



(12) 
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The first condition ensures that transport between scattering events may be treated semiclassically. This limit can 
be achieved for our system when k — > oo, T — > while the classical parameter if = feT remains constant. When 
ijif < L as long as the motion is not localized (i.e. L <C £,) it is diffusive, since a particle scatters many times before 
it can traverse the system. The resulting mean free path for our system (QJ) is Im — ^/Dk (see Appendix I) while the 
localization length for d = 1 is £ ~ D k /2 6] and for d = 2 is £ ~ l M e Dk/ * Hll^- 

A great advantage of the kicked rotor consists in the fact that due to the unitarity of the evolution operator all 
its eigenstates have the same statistical properties. This is in contrast to the eigenstates of Hamiltonian models 
(like Anderson model), where the eigenstates belonging to different parts of the spectrum have different statistical 
properties and therefore they must be picked up from a narrow energy window. This allow us to collect a huge set of 
numerical data and perform a rather accurate numerical analysis. The results of our investigation will be discussed 
in section V. 



IV. THE SCATTERING APPROACH 



To proceed with the analysis of the resonance widths and delay times we turn the closed 2D KR model (0 into an 
open one. To this end we impose 'absorption' at the boundary of a square sample of size L x L in the momentum 
space in complete analogy with the classical set-up. In other words, every time that one of the components of the two 
dimensional momentum [Ci, £2) takes on the value 1 or L, the particle leaves the sample. The corresponding unitary 
scattering matrix S (see Appendix III) is given [24[ 

S(lo) = a/7 - WW* - WU 1 (13) 

- VI - W^WU 

where I is the L 2 x L 2 unit matrix and W is a M x L 2 matrix with matrix elements 

Wij(k,A)={ ' '' (14) 
0, i?j 

with wf £ [0, 1] the tunnel probability in mode i. In the case of perfect coupling, which is considered here, Wi = 1. 
Then WW^ — Imxm and W^W is a L 2 x L 2 diagonal matrix with M non-zero elements equal to one. From the 
physical point of view, W describe at which "site" of the Lx L sample we attach M "leads" [in our case M = 4(L — 1)]. 
Here W^W is a projection operator onto the boundary, while P = I — W^W is the complementary projection operator 
which satisfies 



\Jl- W^W = I -W f W (15) 
Taking this into account, the expression (|13f) can be simplified: 

S(lu) = -WU — : — ; — —W^ = WU — j W^. (16) 

and using the unitarity of the evolution operator U we can rewrite the last expression as follows: 

S = W (17) 

The scattering matrix (|17|l can cquivalently be written in the form used conventionally in quantum chaotic scattering: 

S(cj) = -WUe iu - 1 — 7T W^ (18) 

The scattering matrix Sij given by Eq. (|18(l can be interpreted in the following way: once a wave enters the sample, 
it undergoes multiple scattering induced by [I — e^PU]^ 1 — Y^=<d ( e *" PU) n until it is transmitted out. It is clear 
therefore that the matrix U = PU propagates the wave inside the sample. However, contrary to the closed system in 
which the evolution operator is unitary, the absorption breaks the unitarity of the evolution matrix U. 

The eigenvalues of U occurring at complex quasi-energies u> n = uj n — ^T n are the poles of the scattering matrix. 
They represent long-lived intermediate states to which bound states of a closed system are converted due to coupling 
to continua. Here Lo n and T n are the (dimensionless) position and width of the resonances, respectively. 
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Having at our disposal the scattering matrix S we can calculate the Wigner-Smith delay time. It captures the 
time-dependent aspects of quantum scattering and formally it is defined as 

tw{uj) = -i— bidet S(u;) (19) 
dui 

and can be interpreted as the typical time an almost monochromatic wave packet remains in the interaction region. 

A generalization of the notion of Wigner delay time is given through the Wigner-Smith operator. The latter is 
defined as (for the kicked rotor model one should use quasi-energy instead of energy in the definition of Wigner-Smith 
operator) 

Q = -S^ (20) 
i ouj 



Introducing a new operator 



K= — t j — - (21) 



and taking the derivative of both sides of Eq. I|17(l we obtain 



= -ie-^WKU^KW^ (22) 



Then the definition of the Wigner-Smith yields 

Q = -ie-^WK^WWKU^KW^ (23) 

The Wigner delay time can be expressed as the sum of proper delay times T q . The latter are the eigenvalues of the 
Wigner-Smith operator (|23[) . 



V. DISTRIBUTION OF EIGENFUNCTION INTENSITIES 



The statistical properties of wavefunction intensities have sparked a great deal of research activity in recent years. 
These studies are not only relevant for mesoscopic physics El El El El El IT^JTl El IU |2a Eil |H El H3 , 
but also for understanding phenomena in areas of physics, ranging from nuclear [31| and atomic |32tl33| to microwave 
physics |3 HI HI Ell an d optics [38| . Experimentally, using microwave cavity technics it is possible to probe the 
microscopic structure of electromagnetic wave amplitudes in chaotic or disordered cavities j3J, |2a, l36l l37| . Recently, 
the interest in this problem was renewed when new effective field theoretical techniques were developed for the 
study of the distribution of eigenfunction intensities 7- > (|?/>| 2 ) of random Hamiltonians. As the disorder increases, 
these results predict that, the eigenfunctions become increasingly non-uniform, leading to an enhanced probability of 
finding anomalously large eigenfunction intensities in comparison with the random matrix theory prediction. Thus, 
the notion of pre-localized states, i.e. states which are localized much more strongly than typical eigenstates, has been 
introduced 

El El El El El E3 to explain the appearance of long tails in the distributions of the conductance and 
other physical observables [lj ■ 

Up to now all theoretical predictions ElElElElElElElEHHEHHand numerical calculations [HIHHI 
apply to disordered systems and are based on an ensemble averaging over disorder realizations. Their validity, however, 
for a quantum dynamical system (with a well defined classical limit) that behaves diffusively is not evident. Here we 
show that pre-localized states exist also for dynamical systems with underlying classical diffusion and investigate the 
effect of these states in "P(|t/>| 2 ). An example of them is reported in Fig. (0J. Our main conclusion is that in a generic 
dynamical system with classical diffusion, 7 , (|'0| 2 ) is described quite well by the nonlinear a— model (NLSM). We 
point out that between the various theoretical works there is a considerable disagreement about the parameters that 
control the shape of 7- > (|V'| 2 ) and their dependence on time-reversal symmetry (TRS). More specifically, the NLSM 
suggests that the tail of 'P(|'0| 2 ) in two dimensions (2c?) is sensitive to TRS [la. fl7i lla . fl9l l25|| . while a direct optimal 
fluctuation (DOF) method predicts a symmetry independent result |20| . This prediction was an additional motivation 
for the present study. 

In the next subsection we calculate the distribution function V(y = L d \^k{n)\ 2 ) by using a direct diagonalization 
of the Floquet operator ©. The TRS is broken entirely for a = 5.749, i.e. the phase difference for a typical trajectory 
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FIG. 4: An example of anomalously localized state. The size of the system is L x L = 90 x 90, the diffusion coefficient D = 33.8 

and its time-reversal counterpart is larger than 27r, so that all interference effects between time-reversal trajectories 
are suppressed. In order to test the issue of dynamical correlations, we randomize the phases of the kinetic term of 
the evolution operator © and calculate the resulting V(y). This model will be referred to as Random Phase KR 
(RPKR). Since all our cigcnfunctions have the same statistical properties (in contrast to the Anderson cases where 
one should pick up only eigenfunctions having eigenenergies within a small energy interval |27ll28ll29| ) we make use 
of all of them in our statistical analysis. The classical parameter K > 6.36 is large enou gh i n all cases to exclude the 
existence of any stability islands in phase space, the accelerator modes known for KR |2ll l3^ | are avoided as well. 
The influence of the accelerator modes on the eigenfunction statistics is an interesting problem, which is out of the 
scope of the present study and deserves a separate investigation. The classical diffusion coefficient Dk is calculated 
numerically by iterating the classical map obtained from IH. Below we present our numerical results and compare 
them to the predictions of Refs. E3, Q EH E3, E3, EM. E3. E3, - 



A. Numerical results 



Let us start our numerical analysis from the ballistic regime where g — > oo. In this case, RMT is applicable and 
one finds 0: 

V^(y) = cxp(-y/2)/V2^ (24) 

-P(p=2)(y) = exp(-y). (25) 

Here (3 denotes the corresponding Dyson ensemble: (3 — 1(2) for preserved (broken) TRS. This result can be easily 
understood. Indeed, within the random matrix theory one assumes that all the eigenvector components are inde- 
pendent (the normalization of the eigenvector is not essential in the thermodynamic limit, i.e. when the number of 
its components becomes very large) random variables obeying Gaussian distribution. Going to the distribution of 
the modulus square of the components one immediately recovers Eq. (|24() . For the case of the broken time reversal 
symmetry one should take into account that each component has statistically independent real and imaginary parts, 
leading to the distribution given by Eq. (j25(l . The numerical data presented in Fig. © shows that the distributions 
of the eigenfunction intensities in the ballistic regime for the two-dimensional kicked rotor are described very nice by 
the RMT prediction. 

As the ratio between the system size and the localization length increases, the deviations from the RMT results of 
the body and the tails of the distribution Vpiy) become noticeable and can be parameterized by a single parameter 
which is the dimensionless conductance g = Dk- For y < Dk, according to all studies p^4L fl5L [la. fl7L Il8f| V(y) is just 
the RMT result with polynomial corrections in powers of L/Dk, i.e. Vpiy) = 'Pg' T (y)[l + 8Vp(y)]. The leading term 
of this expansion is given by 



(26) 
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FIG. 5: Distributions of the eigenfunction intensities in the ballistic regime for two-dimensional kicked rotor (solid lines) 
compared with the RMT predictions (dashed lines). 



where n is the 2d diffusion propagator (time- integrated return probability), which is identical for (3 — 1 and (3 = 2 
since it is a classical quantity. 

Figures EJ a ) and (b) show corrections to Va' MT for g = Dk 3> 1 for two representative values of Dk- We find again 
that the form of the deviations are very well described by Eq. (|26|l and the agreement becomes better for larger values 
of the diffusion constant. This is due to the fact that by increasing Dk we are approaching the semiclassical region 
and therefore Eqs. (|12[1 are better satisfied. At the same time higher order corrections in SVfj(y) become negligible 
with respect to the leading term given by Eq. (|26[) . 

In Fig.[Sfc) we summarize our results for various Dk values. The extracted up values are obtained by the best fit of 
the data to Eq. 1)26(1. Again we find that up depends linearly on 1/Dk- However, n\ and K2, are different. Moreover 
the best fit with n fj = ApD^ 1 + Bp yields Ap=i = 5.44 ± 0.03 and A fi=2 = 10.84 ± 0.04 indicating that the ratio 
R = A 2 /Ai is close to 2, a value that could be explained on the basis of ballistic effects 0, |25l |53, . Taking the 
latter into account leads to an additional term in the classical propagator up — Kdiff + ^^baii- The first term is the 
one discussed previously and is associated with long trajectories which are of diffusive nature while the latter one is 
associated with short ballistic trajectories which are self-tracing 0,|2j|. Thus, when Kdiff "C Kbaii we get R = 2. 
The calculation with the RPKR model shows, however, that the corresponding ratio is R ~ 1 in agreement with the 
theoretical prediction for disordered systems with a pure diffusion. This indicates that dynamical correlations can be 
important. 

For the tails of the distributions, the result of the NLSM within a saddle-point approximation 0, 0, 0, |2^| is 

^(y^exphC^lny) 2 ], g|= ^ 2p ^ (27) 

Note that the decay in the tails of Eq. (|27|l depends on (3. Recently, a DOF method was used to calculate the tails 
of Vp(y) |20| for the white- noise random potential. It was found that the tails are still given by Eq. I|27|) but with a 
log-normal coefficient C which is independent of the parameter (3 : 

c =1Tp HLf\y (28) 

Figure 0a) shows a representative case of Vp^i^y > Dk). The tails show a log-normal behavior predicted by 
Eq. H27fl • In Fig. 0b) we report the log- normal coefficients Cp extracted from the best fit to our numerical data, 
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FIG. 6: Corrections to the distribution intensities 5Vp{y) for the two- dimensional kicked rotator model ®. The system size 
is L = 90, (o) p = 1, (o) f3 = 2. The solid (dashed) lines are the best fit of J^HJ for /3 = 1(2) to the numerical data: (a) D k ~ 34 
and (b) Dk ~ 53 ; (c) Fit parameters up vs. D^ 1 . The solid (dashed) lines are the best fits to Kp = ApD^ 1 + Bp for (3 = 1(2). 



versus the classical diffusion coefficient. A pronounced linear behavior is observed in agreement with both theories. 
However one clearly sees that C\ differs from Co in contrast to the DOF prediction (|28|) and to recent numerical 
calculations done for the 2c? Anderson model (23. We point out here that in j2^| the authors were not able to go 
to large enough values of conductance g (in comparison to our study) where the theory can really be tested. In 
contrast, the NLSM predicts a value of 2 for the ratio R = /Cf . We note that Cp is only the leading term in 
Dk- In order to calculate this ratio, we performed a fit to our data with Cp = ApDk + Bp. The resulting ratio was 
found to be R — A2/A1 — 1.97 ± 0.03 in perfect agreement with the NLSM predictions. Finally in Fig.[7Jb) we also 
present our results for the RPKR model (using the same data as the one in Fig. Etc)). Again we found that the ratio 
R = 1.96 ±0.03 rj 2. Thus V(y > D k ) depends on TRS and is described by the NLSM. The fact that the prediction of 
a DOF method is not observed in our calculations might be due to non-universal (depending on the type of disorder) 
character of this result. 



VI. OPEN SYSTEMS: DISTRIBUTION OF RESONANCES AND DELAY TIMES 

In this section we analyze the statistical properties of resonances V(T) and Wigner delay times V(t) for a ID 
diffusive system. This study is important for various applications among which are diffusive random lasers and 
microwave cavities 0, ^2 Hd> where most of the theoretical treatment is limited by RMT. We point here that 
current developments of microwave experiments in random dielectric media in the diffusive regime 45] may allow us 
a direct comparison between theory and experiment. 

Specifically, we have found that the resonance width distribution 'P(r) is given by 

V(T<r c i) - exp(-C /3 (lnr) 2 ) , where C fj ~f3D 

PFZrj-M* (29) 
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FIG. 7: (a) Tails of the distribution Vp=\{y > D k ) for the model Q and D k ~ 35. The system size is L = 80, (o) /3 = 1, (o) 
j3 = 2. The solid (dashed) lines are the best fit of 12711 for (5 — 1(2) to the numerical data; (b) Fitted log-normal coefficients 
C/3 versus the classical diffusion coefficient Dk- The solid (dashed) lines are the best fits to Cp = ApDk + Bp for j3 = 1(2). 

while the distribution of the Wigner delay times, is given by the following expressions: 

-Pir^r-, 1 ) ~ -Lexpi-a/r) 

T(t>T- 1 1 ) ~ exp^C^lnr) 2 ) (30) 

where D is the classical diffusion constant, /? denotes the symmetry class and a is some constant of order unity. 

Our theoretical considerations are supported by numerical calculations for the open analog of the 2D KR model 
described in sections II and IV. The parameters of the model were chosen in such a way that the conditions H12|) 
discussed in the previous section were fulfilled. In order to improve our statistics, we randomized the phases of the 
kinetic term of the evolution operator and used a number of different realizations. The results of the previous 
section allow us to conclude that this procedure doesn't change the universal features of the model. In all cases we 
had at least 60000 data for statistical processing. 



A. Resonance widths distribution 

For diffusive mesoscopic samples, there is no systematic investigation of V(T) besides Ref. 0] where the authors 
have focused on the large tails of V{T) for a quasi-lD system in the diffusive regime. This deficiency is felt especially 
strong in the random-laser community where one wants to know the statistical properties of the lasing threshold. 

Usually we model random laser as a disordered or chaotic system containing a dye that is able to amplify the 
radiation with a rate rj, in a certain frequency interval. In contrast to the traditional lasers where the necessary 
feedback is due to mirrors at the boundaries of the laser cavity, the key mechanism for random lasers is the multiple 
scattering inside the medium |4£j. The lasing threshold is give n by the value of the smallest decay rate (i.e. smallest 
resonance width) of all eigenmodes in the amplification window [47ll48j| . The underlying reasoning is that in the mode 
with the smallest decay rate the photons are created faster by amplification than they can leave (decay) the sample. 
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Assuming that the number of modes K ^> 1 that lies in the frequency window where the amplification is possible, 
have resonance widths T that are statistically independent one gets for distribution of lasing thresholds V(T) pTllisl 
IH: 



V(T) = KV(T « 1) 



i - / V(V « i)dr' 

Jo 



K-l 



(31) 



where we have assumed that all K resonances are distributed according to ^(r << 1). The validity of this approxi- 
mation was verified recently in the framework of the RMT 50] . An important outcome of our study will be that one 
can identify traces of pre-localized states in the latter distribution and consequently in V(T). This send some light to 
recent experimental finding for random lasers which suggests the appearance of localized modes in diffusive samples 

We start our analysis with the study of resonance width distribution V(T) for T < T c i. The small resonances T < A 
can be associated, with the existence of pre-localized states of the closed system which were discussed in the previous 
section. They consist of a short-scale bump (where most of the norm is concentrated) and they decay rapidly in a 
power law fashion from the center of localization [20l |25| . One then expects that states of this type with localization 
centers at the bulk of the sample are affected very weakly by the opening of the system at the boundaries. In first 
order perturbation theory, considering the opening as a small perturbation we obtain 

l = {*\w*wm= J2 i*wi 2 ~£i*wi 2 (32) 

nGboundary 

where |^(L)| 2 is the wavefunction intensity of a pre-localized state at the boundary. At the same time the distribution 
of 9 = l/\/L^(L) for large values of the argument is found to be of log-normal type pfij ]: 

V{9) ~ exp (-tt 2 D \n(6 2 )) (33) 

Using this together with Eq. I|32(l we obtain 

T(l/T) ~exp(-7r 2 Dln 2 (l/r)) (34) 

We would like to stress that, based on the results of the previous section, the expression for V(6), must be corrected 
by including the TRS factor (3 in the exponent. Taking all the above into account we end up with the expression 
given in Eq. J2HJ). 

The numerical data reported in Fig. |H1 support the validity of the above considerations. However, we would like to 
mention that the perturbative argument is valid only for the case of very small resonances i.e. T < A, whereas our 
numerical data indicate that one can extend the log-normal behavior of V(T) up to resonances with A < T < T c /. 

Next we turn to the analysis of V(T) for T > r c /. In Fig.^a) we report our numerical results for V(T) with preserved 
(broken) TRS for two representative values of D. An inverse power law V(T) ~ r -1 ' 5 is evident in accordance with 
Eq. I|29|) . (The behavior of the extreme large T tails of V{T) is essentially determined by the coupling to the leads 
which is model dependent. Their relative number is proportional to M/ L 2 ~ L^ 1 and therefore they are statistically 
insignificant.) ^From the figure it is clear that this part of the distribution is independent of the symmetry class, in 
contrast to the small resonance distribution discussed above. 

The following argument provides some understanding of the behavior of V(T) for T > r c /. First we need to recall 
that the inverse of T represents the quantum lifetime of a particle in the corresponding resonant state escaping into 
the leads. Moreover we assume that the particles are uniformly distributed inside the sample and diffuse until they 
reach the boundaries, where they are absorbed. Then we can associate the corresponding lifetimes with the time 
tn ~ l/r^j ~ R 2 /D a particle needs to reach the boundaries, when starting a distance R away. This classical picture 
can be justified for all states with F > T c i ~ D/L 2 . The relative number of states that require a time t < in order 
to reach the boundaries (or equivalently the number of states with V > Tr) is 



/>oo 

Pmt(r«) = / v(r)dr 

JTr 



s(t R ) 



L 2 

where S(t]i) is the area populated by all particles with lifetimes t < t^. In the case of open boundaries we get 



(35) 
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o 



ln 2 (l/r) 




FIG. 8: (a) The distribution of resonance widths (plotted as V(l/T) vs. 1/F) for T < F c i for two representative values of D. 
The system size in all cases is L = 80. Filled symbols correspond to broken TRS. The solid lines are the best fit of Eq. 129^ for 
j3 = 1(2) to the numerical data, (b) Coefficients Cp vs. D. The solid lines are the best fits to Cp = ApD + Bp for f3 = 1(2). 
The ratio R = AijAx = 1.95 ± 0.03 



For Tr > Yd the first term in the above equation is the dominant one and thus Eq. 1)29(1 follows. 

Here it is interesting to point that a different way of opening the system might lead to a different power law behavior 
for 'P(r). Such a situation can be realized if instead of opening the system at the boundaries we introduce "one-site" 
absorber (or one "lead") somewhere in the sample. In such a case we have 



V int (T R ) 



S(t R ) 
L 2 



L 2 



DtR 

L 2 



rV 



(37) 



The above result is valid for any number M of "leads" such that the ratio M/L 2 scales as 1/L 2 . In Fig. EJb) we 
report the integrated resonance width distribution 'Pint(r) for the case with nine "leads" attached somewhere to the 
2D sample. 

A straightforward generalization of our arguments for 3D systems in the diffusive regime gives 

/r^ 

Jt r T r + 3 ^ 



V int (T R ) ~ V - 2 7^ + 3 [j?) (38) 
which for T R > T c i leads to the same universal expression as in Eq. (|29|l . Similarly, the analog of Eq. (|37|l in 3D is 



' ( ^ 



(39) 



It is interesting to compare the above prediction 1(29(1 with the results of the random matrix theory. In the general 
case, Fyodorov and Sommers [52^ proved that the distribution of scaled resonance widths 7 = T/A for the unitary 
random matrix ensemble, is given by 



yi> r(M ) 7 



(h. 



_ 7wg smh(77r)\ 

( 7 7T) J 



(40) 
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FIG. 9: (a) The resonance width distribution V(T) for preserved TRS and D = 20.3 (o) and D = 33.5 (o). The corresponding 
filled symbols represent V(T) for broken TRS and the same values of D. The dashed (solid) vertical line mark the classical 
decay rate r c j for D — 20.3(D = 33.5). (b) The Vi n t{T) for a sample with nine leads (lower curve). For comparison we plot 
also the Vint(T) for the same sample but when we open the system from the boundaries. The dashed lines correspond to the 
theoretical predictions 136H and 1371 . 



where M is the number of open channels and the parameter q controls the degree of coupling with the channels (q = 1 
for perfect coupling). Our numerical data from the 2D KR model in the ballistic regime, reported in Fig. 1101 are in 
excellent agreement with the theoretical prediction H40(l . 

In the limit of M ^> 1, which is relevant for the comparison with Eq. (|29[) . Eq. (|40|l reduces to the following 
expression 52] 

f _M_ f or m < ■y < M 

U} \ otherwise ' V ; 

One can see two essential distinctions between this result and Eq. (|29[l . Firstly, the power law P(T) ~ 1/r 2 is not the 
same as the power law predicted by Eq. I|29[l for large resonances P(T) ~ 1/r 3 / 2 . Since this difference appears in the 
"classical" part of the distribution, it can be explained as a difference in the classical dynamics of a particle inside 
the system: ballistic (RMT) versus diffusive motion. Indeed, taking into account that for ballistic system R ~ vt and 
r c ; ~ v/L, where v is the velocity of the particle, one immediately finds from Eq. 136|) that VintO^R.) ~ ^ci/^R for 
r.R > r c /, in agreement with the RMT prediction P(T) ~ 1/r 2 . Secondly, according to Eq. (|4lT) there is a gap in the 
distribution of the resonance widths: there are no resonances with widths smaller than , M , . , . The existence of the 
gap can be understood, if one relates the small resonances to the coupling of the wavefunctions to the leads. Since 
the wavefunctions in the RMT are extended, the probability to find a wavefunction, which is weakly coupled to all M 
channels, goes to zero when the number of channels becomes very large M ^> 1. In the diffusive regime, in contrast, 
there are pre-localized states, which are weakly coupled to the leads. Due to their existence the distribution of the 
small resonance widths has a non-trivial behavior described by Eq. I|29|) . 
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FIG. 10: Distributions of the resonance widths in the ballistic regime for two-dimensional kicked rotor (o) compared with the 
RMT predictions (dashed lines). 



B. Wigner delay times distribution 



We turn now to the analysis of Wigner delay times and calculate their probability distribution. It can be shown 
that this distribution is related to the distribution of reflection coefficients R in the present of weak absorption. 
Absorption is one of the main ingredients in actual experimental situations and its theoretical understanding is of 
great importance. Unfortunately a comprehensive treatment of absorption is still lacking. There are only very few 
reported analytical results for the distribution of the reflection coefficient V(R) in the presence of absorption and all 
of them are either within the regime of applicability of RMT EH El El El Hi or applied for quasi- ID geometry 
|53l l54| . Here we will derive V{R) in the diffusive regime in two dimensions and in the weak absorption limit. In this 
limit it was shown El EH that the following relation between the proper Wigner delay times and reflection coefficients 
holds: 

R q = 1 - T q /r a (42) 

where r q are the proper delay times (eigenvalues of the Wigner-Smith operator) and l/r a is the absorption rate. Thus 
the knowledge of 'P(R) reduces to the calculation of the distribution of proper Wigner delay times V(j q ) |42l |. 

Bellow we make the standard assumption that the resulting distribution generated over different energies is equiv- 
alent with the one generated over different disorder realizations. Our starting point is the well known relation 

L 2 

T H = E 7 ^Tr2 77 ( 43 ) 

which connects the Wigner delay times and the poles of the S— matrix. 

Let us start with the far tails. It is evident that large times t{oj) ~ T^ 1 corresponds to the cases when ui ~ u> n 
and r„ <C 1. In the neighborhood of these points, t(lo) can be approximated by a single Lorentzian l|43|l . Sampling 
the quasi-energies u> with step Auj <§; T m i n we calculate the number of points for which the time delay is larger than 
some fixed value r (see Fig. II 1|) . Assuming that the contribution of each Lorentzian is proportional to its width one 
can estimate this number as J2r <i/r r n /Ai£. For the integrated distribution of delay times in the limit Auj — > we 
obtain 

,1/t 

V int {r) ~ / dTV(T)T (44) 
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FIG. 11: Schematic plot for the Wigner delay time as a function of quasi-energy according to Eq. I|48|l . 

and by substituting the small resonance width asymptotic given by Eq. I|29[) we come out with the log-normal law of 
Eq. H3()|l in agreement with our numerical findings reported in Fig. 1121 

Now we estimate the behavior of V{t) for r < rjj 1 . In this regime many short-living resonances contribute to the 
sum l|43|) . We may therefore consider r as a sum of many independent positive random variables each of the type 
T n = r„j:„, where x n — Suj~ 2 . Assuming further that 5u> n are uniformly distributed random numbers we find that 
the distribution V(x n ) has the asymptotic power law behavior 1/xf/ 2 . As a next step we find that the distribution 

3/2 3/2 

Vijn) decays asymptotically as l/r n where we use that V{T n ) ~ ■ Then the corresponding V(t) is known 

to be a stable asymmetric Levy distribution L m ,i(t) of index [i = 1/2 55] which has the form given in Eq. I|3(JII at 
the origin. We point out here that the asymptotic behavior V{r) ~ 1/r 3 / 2 emerges also for chaotic /ballistic systems 
where the assumption of uniformly distributed Suj n is the only crucial ingredient (see for example |52(). 

Since t = X^=i T qi we expect the behavior of the distribution of proper delay times V[r q ) to be similar to V{r) for 
large values of the arguments (for t>1 we have r ~ r™ ax ). Moreover, from the numerical point of view V(r q ) can 
be studied in a better way because a larger set of data can be generated easily. Our numerical findings for V[r q ) are 
reported in Fig. 1121 and are in nice agreement with Eq. even for moderate values of T q . We stress here that the 
dashed lines in Fig. 1121 have slopes equal to Cp taken from the corresponding log-normal tails of V(T). 

Finally we would like to compare our result 13U[) with the results known from the random matrix theory. Although 
the distributions of the proper delay times 15a. l57| and partial delay times (defined as a derivative of the partial phase 
6i of the 5*-matrix with respect to energy) 52] have been calculated recently, there is no explicit analytical formula for 
the distribution of the Wigner delay times for number of open channels M > 2 [5a . Nevertheless using consideration 
similar to one presented above (see the discussion for r < T^ 1 ) one can argue [52] that the part oIV{t) for r < T^ 1 is 
the same in the RMT as stated in Eq. (|3L)fl . However, the distribution of the large delay times in the RMT is expected 
to have a power law behavior V{r) ~ l/ r 2 +/ 3M / 2 w ith M being the number of open channels. This is in contrast with 
a log-normal tail stated in Eq. (f 301) . 

VII. CONCLUSIONS 

This paper is devoted to the investigation of chaotic and disordered systems characterized by the classical diffusion. 
Section V deals with study of the closed system. Specifically, we perform a detailed numerical analysis of the eigen- 
function intensities V(y) of the two-dimensional kicked rotor on a torus. Our results indicated that the distribution 
V{y) of generic quantum dynamical systems with diffusive classical limit is affected by the existence of pre-localized 
states. The deviations from RMT are well described by field theoretical methods developed for disordered systems. 
In particular, we find that the dependence of the tails of Vp(y) on TRS is described correctly by NLSM. 

Section VI deals with the study of the open system. After introducing a scattering formalism for the KR model 
we investigated the distribution of the resonance widths V(T) and Wigner delay times V(t). We obtain the forms of 
these distributions (log- normal for large r and small T, and power law in the opposite case) for different symmetry 
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FIG. 12: The proper delay times distribution V{j q ) for D = 20.3 (o) and D = 29.8 (□ ). The (•) correspond to D = 20.3 but 
now with broken TRS. The dashed lines have slopes equal to Gp extracted from the corresponding V(T) (see Fig. 8). In the 
inset we report V(r q ) for moderate values of r 9 in a double logarithmic scale. 

classes and show that they are determined by the underlying diffusive classical dynamics and by the existence of the 
pre-localized states. Our theoretical arguments are supported by extensive numerical calculations. 

Although the arguments, we used to explain the behavior of T(T) and V(t), can be easy generalized to the three- 
dimensional case, the numerical test of these predictions has not been still performed. Moreover the study of three- 
dimensional case would allow to investigate these distribution at the critical point of the metal-insulator transition. 
The first attempt to attack this problem was done in Ref . |59j| . but a detailed understanding is still required. 

We thank I. E. Smolyarenko for useful comments. This work was supported by a grant from German-Israeli 
Foundation for Scientific Research and Development. 

VIII. APPENDIX I: THE MEAN FREE PATH 

Let us consider the evolution operator of the two-dimensional kicked rotor introduced in Eq. © (to simplify the 
calculations we put T = 1) 

U _ e -iH (CuC2)/2 e ~iV(e 1 ,e 2 ) e -iH (c 1 ,C2)/2_ (45) 
The set of the orthogonal eigenfunctions of the free Hamiltonian Hq is given by plane waves in the angle representation 

4>(n 1 ,n 2 )(0i,e 2 ) = l-e i ^ +n ^. (46) 

They are normalized in such a way that 

p2tT flit 

/ / dd 1 de 2 \c/ )(nun2) (9 1 ,e 2 )\ 2 = 1. (47) 

Jo Jo 
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Let us denote by Un,^^\< ni ,n 2 ) the matrix elements of the evolution operator in this basis 

U (ki,k 2 ),(ni,n 2 ) = <0(fc 1 ,fc 2 )|C / l0(n 1 ,n 2 ))- 



(48) 



According to the definition of the evolution operator the modulus square of its elements have the meaning of the 
probability to change the initial momentum (nj,n 2 ) to the final momentum (ki,k 2 ) in one kick. Therefore one can 
define the mean free path in momentum space |8| by 



l M = E +r 2) 1^(0,0), (n.rj 



(49) 



Here we used the fact that \U ikl+miM+m2)Ani+mun2+m2) \ 2 = |^(fc 1 ,fc 2 ),(«i,n 2 )| 2 , so without loss of generality one can 
take the initial momentum equal to (0, 0). In order to calculate the right hand side of Eq. I|49(l we first give an explicit 
expression to the matrix elements C^(fe 1 ,fc 2 ),(ni.n 2 ) : 



(fei,fc 2 ),(ni,rs 2 ) 



i,j,s,t 



-i(H (k 1 .k 2 )+H Q (n 1 .n 2 ))/2 _j_ f * dQ^Q^-Wtfxfa) e i(m-ki)9i+i(n 2 -k 2 )9 2 



Taking (k\,k 2 ) = (0,0) and (nx,n 2 ) = (ri,r 2 ) we obtain 



U, 



(o,o), ( 



= e -i(H ( ,0)+H (r u r 2 ))/2_L f*" f*" ^Q^V (9 U 9 2 ) ^9, jr 2 9 2 
47T Jo Jo 



The substitution of this expression into Eq. I|49|) yields 



(50) 



(51) 



2ir r 2% r 2n r 2n 



o Jo Jo 



Taking into account that 



,2ir l (9 1 -91 l ) _ 



If 



iri(9 x -9lx) 



d8 1 d9 1 

and that the same is valid for r 2 the partial integration of the Eq. l|52|l gives 



(53) 



;2 



(4tt 2 ) 2 



27T /.27T r 2n i-2-K 



JO JO JO 



d0 x d0 2 d0 x d0 2 



2 



e -i(V(9 u 9 2 )-V{ei,9 2 )) e iri(e 1 -e 1 ) e ir2(9 2 -9 2 ) 

The summation of the exponents over r\ and r 2 yields two (^-functions: 

r 2TT n 2ir ?2ix r 2ir / q2 f)2 



M 



1 

4^2 



o Jo Jo Jo 



d0 1 d0 2 d0 1 d0 2 



t) 2 



hdOi de 2 de 7 



e 

I 

4^2 



i{v(9 U e 2 )-v(9 U h)) 5 (e x _ fygfo _ § 2 ) = 

2ix /-27T r 2n r 2n / dV 



d 2 



dQxdQi 80 2 d0 2 



o Jo Jo Jo 

i(V(9 1 fi 2 )-V(9 1 ,9 2 )) 



de 1 de 2 de 1 de 2 [ ^--(e 1 ,d 2 )^-(§ 1 ,e 2 ) + ^-(d 1 ,e 2 )^-(e 1) e 2 ) 

ooi d6\ ov2 o0 2 



6(0! ~ 1 )S(0 2 - 2 ) = 



1 

47T 2 



d6id6 2 



dV \ ( 8V 



(54) 



(55) 
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The last expression can be written in a compact form 



I 2 - J- 



Now we calculate the mean free path in the case where the potential V(9i, 62) is given by Eq. iffll: 



V(0i, 6> 2 ) = fc I cos(^) cos(6> 2 ) cos(a) + - sin(20i) cos(26» 2 ) sin(a) 

Taking the derivative of this expression with respect to 6\ and 62 one has 

2 



<9U 

ST*"* 



au 
ae 2 



(0i 



(56) 



(57) 



= A; (sin 6>i cos 9 2 cos a + cos 29 1 cos 26> 2 sin a - 
— 2 sin 9i cos 2#i cos 62 cos 2# 2 cos a sin a) 

2 

= A: 2 (cos 2 #1 sin 2 6> 2 cos 2 a + sin 2 26 1 sin 2 26* 2 sin 2 a + 

+2 cos 9\ sin 2#i sin 6* 2 sin 2# 2 cos a sin a) (58) 

The integration over 9\ and # 2 yields for the mean free path 

1 k 2 
l 2 M = -^k 2 (7r 2 cos 2 a + 7r 2 sin 2 a + tt 2 cos 2 a + tt 2 sin 2 a) = — (59) 

which is the same result with the one derived in || for different potential V (61,62) ■ 

IX. APPENDIX II: DIFFUSION COEFFICIENT IN THE RANDOM PHASE APPROXIMATION 



Here we give for completeness of the presentation the derivation of the diffusion coefficient for our model (11121) in 
the random phase approximation. We start by writing the classical maps Eq. for the general form of the potential 
V(6 X ,6 2 ): 

6i(n+l) = 6i(n) + TiTCx(n) mod 2tt 
6 2 (n + l) = 6 2 (n) + T2TC 2 (n) mod 2tt 

CAn+1) = C 1 (n)-^-(6 1 (n + l),6 2 (n+l)) 

dV 

£ 2 (n + l) = C 2 (n)- — (6i(n + l),6 2 (n + l)) 

Ot>2 

~ sin(20i(n + 1)) sin(26» 2 (n + 1)) sin(a)) 
The diffusion coefficient in momentum space is defined as 

(C 2 (t) + C 2 (t)) 



D = lim 

t— >oo 



t 



(60) 



(61) 



The average in this expression is taken over an ensemble of trajectories, with different initial conditions. Using the 
classical maps for L\ and £ 2 we obtain 



D = lim — 

n^oo nT 



£ ( (f^ i( ^ 2(i)) ) 2 + {w 2 m)Mi)) 



n n 



E E (^(»iWA(i))^(fli(?),«2C7)) + ^(<5'iW > fl2(i))^(fliO'),fe(7))) 

i=l j=lj^i 



(62) 



For the large values of the kicking strength k, in a good approximation one can consider the phases 6\(i) and 6 2 (i) 
as random variables which are uncorrelated for different i and distributed uniformly in the interval [0, 2tt]. Using this 
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random phase approximation it is easily to show that only the diagonal terms in Eq. (|62() give non-zero contribution 
in the limit n — > oo. Taking into account that the distribution of the phases is uniform one can convert the sum for 
the diagonal terms into an integral. Finally we obtain 




Then the formula for the diffusion coefficient measured in number of kicks Dk — TD is given by 

Dk = i C fj ddM ^ ¥{du 02)1,2 (64) 

which has exactly the same form as one appearing Eq. i(56() . Thus we obtain that in the random phase approximation 
the following relation between the mean free path and the diffusion coefficient is valid: 

Dk = 4 = y ( 65 ) 

Therefore changing the kicking strength we can easy tune the diffusion constant or Thoulcss conductance (for disor- 
dered systems). This allows us to investigate various regimes: ballistic, diffusive and localized. 



X. APPENDIX III: UNITARITY OF THE S-MATRIX 



Let us rewrite the expression for 5-matrix ((16(1 in a more symmetric way. To this end we use a series expansion for 
inverse operator in 116(1 : 



[U-e-™- W^WU] 1 = 
1 



[U-e 



I- 



1 



U-e- 



-W*WU) 



U 



-W^WU 



1 



U-e- 



E 

k>0 



U-e- 



-W*WU) 



U 



Substituting this expansion in Eq. 1(16(1 we obtain 



S 



i 



U-e- 



-W ] WU ) 



k>0 



J u 



w 



x —^w^ = y( wu — 

k>0 v 7 



fc+i 



U-e- 



W 



U-e 



I-WU V X — Ws 



w[i 



w 



U-e 



Now using the unitarity of the evolution operator U we can calculate the Hermitian conjugate S'-matrix: 



of " e 



W-St 



w- 



rPFt 



w 



U-e 



(66) 



(67) 



(68) 



Thus the unitarity of the 5*-matrix is proven. 
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